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ABSTRACT 

Particulate  composites  consist  of  individual  particles  of  one  material  dispersed  throughout  and  held  together  by  a 
polymer  binder.  The  mechanical  and  physical  properties  of  the  composite  depend  on  the  mechanical  and  physical 
properties  of  the  individual  components;  their  loading  density;  the  shape  and  size  of  the  particles;  the  interfacial 
adhesion;  residual  stresses;  and  matrix  porosity.  Highly-loaded  particulate  composites  are  multi-phase  systems 
that  have  not  typically  been  studied  rigorously,  to  date.  We  are  investigating  whether  or  not  higher-order 
microstructural  features  can  have  a  profound  effect  on  the  static  and  dynamic  mechanical  response  of  these 
multi-phase  (n>2)  polymer-metal-ceramic  composites.  We  present  several  models  for  the  elastic  and  plastic 
behaviors  of  these  materials,  and  compare  the  predictions  with  experimental  data  from  quasi-static  loading 
techniques.  The  high  strain  rate  compressive  properties,  using  a  split  Hopkinson  pressure  bar,  of  these  materials 
are  also  characterized. 

INTRODUCTION 

Particulate  composite  materials  composed  of  one  or  more  varieties  of  particles  in  a  polymer  binder  are  widely 
used  in  military  and  civilian  applications.  They  can  be  tailored  for  desired  mechanical  properties  with  appropriate 
choices  of  materials,  particle  sizes  and  loading  densities.  Several  studies  on  similar  epoxy-based  composites 
have  been  reported  and  have  shown  that  particle  size  [1,  2],  shape  [3],  and  concentration  [4]and  properties  of  the 
constituents  can  affect  the  properties  of  particulate  composites.  In  composites  of  Al203  particles  in  epoxy  (Epon 
828/Z),  increasing  the  particle  concentration  and  decreasing  the  particle  size  were  found  to  increase  the  stress  at 
4%  strain  [5].  A  study  of  aluminum  filled  epoxy  (DGEBA/MTHPA)  found  adding  a  small  amount  of  filler  (~  5  vol.%) 
increased  the  compressive  yield  stress,  but  additional  amounts  of  filler  decreased  the  compressive  yield  stress 
[6].  However,  tests  on  epoxy  (DOW  DER  331/bisphenol-A)  found  that  increasing  the  volume  percent  of  glass 
bead  filler  increased  the  yield  stress  and  fracture  toughness  of  the  material  [7,  8].  In  a  study  on  a  similar  material, 
decreasing  the  aluminum  particle  size  from  micro  to  nano  resulted  in  increased  epoxy  crosslink  density  and 
subsequently  increased  static  and  dynamic  strength  [1].  In  this  study,  a  factorial  design  of  experiments  is  used  to 
examine  the  effect  of  aluminum  particle  size  and  aluminum  and  nickel  volume  fractions  in  epoxy-based  Al-Ni 
particulate  composites. 

EXPERIMENTAL  SET-UP 

Composites  of  aluminum  and  nickel  powders  in  an  epoxy  binder  were  prepared  for  this  study.  The  aluminum 
particle  size  was  varied  between  5  (Valimet,  H5  aluminum)  and  50  (Valimet,  H50  aluminum)  pm.  The  H5 
aluminum,  shown  in  Figure  1(a),  was  found  to  have  an  average  particle  size  of  5.43  pm  with  spherical  smooth 
particles.  The  H50  aluminum  particles,  Figure  1  (b),  were  also  smooth  and  nominally  spherical  with  an  average 
particle  size  of  51.91  pm.  The  nickle  particles  (Micron  Metals)  had  rougher  surface  texture  and  more  irregular 
shape  shown  in  Figure  1  (c),  with  an  average  particle  size  of  47.45  pm.  However,  there  was  also  a  small  fraction 
of  particles  in  the  powder  that  had  an  average  particle  size  of  97.44  pm. 
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Figure  1:  SEM  micrographs  of  (a)  H5  aluminum,  (b)  H50  aluminum  and  (c)  nickel  powders. 

Table  1:  Factorial  design  of  experiments  for  epoxy-based  particulate  composites  containing  aluminum 
and  nickel 


Sample 

Al  Particle  Size 
(pm) 

Al  Volume  Fraction 
(%) 

Ni  Volume  Fraction 
(%) 

MNML-8 

5 

20 

0 

5 

20 

10 

5 

40 

0 

MNML-2 

5 

40 

10 

50 

20 

0 

MNML-3 

50 

20 

10 

MNML-5 

50 

40 

0 

50 

40 

10 

The  composite  materials  used  in  this  study  were  prepared  using  a  fractional  factorial  design  of  experiments,  with 
the  variables  being  the  aluminum  particle  size  (5  or  50  pm),  the  volume  fraction  of  aluminum  (20  or  40  vol.%),  and 
the  addition  of  nickel  (10  vol.%).  The  full  factorial  design  of  experiements  is  presented  in  Table  1,  with  the 
samples  that  were  prepared  highlighted. 

Compression  experiments  at  quasi-static  strain  rates  were  conducted  with  an  MTS  810  testing  system  with  a  100 
KN  test  frame.  Care  was  taken  to  center  the  samples  on  the  platens  prior  to  testing.  MTS  software  was  used  to 
conduct  constant  true  strain  rate  tests  at  6.0  x  10'V1.  A  thin  film  of  Boron  Nitride  (BN)  with  a  layer  of 
Molybdenum  disilicide  (MOSi2)  on  top  was  used  to  lubricate  the  surfaces  of  the  platen  in  contact  with  the  test 
specimen.  Prior  to  testing,  the  compression  frame  went  through  a  warming  up  procedure  of  600  cycles  that 
oscillated  the  hydraulically  driven  ram  at  1  Hz  with  a  displacement  around  2-3  in. 

In  addition  to  the  MTS  system  recording  the  loads  and  displacement  of  the  frame  an  interfacing  software  (VIC 
Gauge  2.0  from  Correlated  Solutions  Inc.)  reads  input  voltages  for  both  the  load  and  displacement.  This  software 
interfaces  with  a  video  system,  which  allows  the  user  to  place  virtual  displacement  gages  on  the  specimen  that 
are  tracked  as  testing  takes  place.  A  high  contrasting  boundary  or  point  is  required  for  tracking,  and  a  black 
marker  was  used  to  draw  fiducial  marks  on  the  specimens.  The  strain  is  measured  directly  from  the  specimen 
itself  rather  than  from  the  MTS  load  frame.  Multiple  virtual  displacement  gages  were  used  for  comparison  and  to 
enable  the  test  to  continue  in  the  event  that  one  gage  failed. 

Compression  experiments  at  intermediate  strain  rates  (103  s'1)  were  conducted  using  a  split  Hopkinson  pressure 
bar  (SHPB)  [9]  system.  The  experiments  were  conducted  using  the  SHPB  system  located  at  AFRL/RWME,  Eglin 
AFB,  FL,  which  is  comprised  of  1524  mm  long,  12.7  mm  diameter  incident  and  transmitted  bars  of  6061-T6 
aluminum.  The  striker  is  610  mm  long  and  made  of  the  same  material  as  the  other  bars.  The  samples,  which 
were  nominally  5  mm  diameter  by  2.5  mm  thick  are  positioned  between  the  incident  and  transmitted  bars.  The 
bar  faces  were  lightly  lubricated  with  grease  to  reduce  friction. 

The  properties  of  the  sample  are  determined  by  measuring  the  incident,  reflected,  and  transmitted  strain  signals, 
si,  £R,  and  eT  respectively,  using  Kulite  AFP-500-90  semiconductor  strain  gages.  These  gages  are  smaller  than 
traditional  foil  gages  and  have  a  much  higher  gage  factor  (140).  The  gages  form  part  of  a  potential  divider  circuit 
with  constant  voltage  excitation,  which  transforms  the  resistance  change  of  the  gages  to  a  voltage  change  and 
compensates  for  temperature  changes.  The  strain  gages  are  dynamically  calibrated  in  situ  by  performing  a 
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number  of  impacts  with  carefully  measured  striker  bar  velocities.  From  the  measured  impact  velocity  and  mass  of 
the  striker,  the  force  amplitude  of  the  stress  pulse  introduced,  F,  can  be  determined  and  compared  to  the  voltage 
output,  V,  from  the  strain  gages  to  give  a  calibration  in  the  form: 


F  =  KV(\  +  bV), 


(3) 


where  K  and  b  are  calibration  factors. 

The  full  derivation  of  the  data  reduction  used  to  calculate  the  strain  rate  and  stress  in  the  specimen,  as  functions 
of  time,  can  be  found  in  reference  [9].  In  order  to  make  representative  measurements  of  material  properties,  it  is 
necessary  that  the  specimen  achieves  mechanical  equilibrium  during  the  experiment,  and  this  is  sometimes 
assumed  as  it  makes  the  strain  rate  calculation  more  straightforward  [9].  The  software  used  in  the  experiments 
presented  in  this  paper  performs  the  one-  and  two-  wave  analyses  automatically  for  every  specimen,  so  stress 
state  equilibrium  is  verified  in  every  experiment.  However,  the  calculation  of  strain  rate  does  not  assume 
mechanical  equilibrium,  rather  it  uses  all  three  of  the  incident,  reflected  and  transmitted  force  pulses  to  calculate 
specimen  strain  rate  through  the  following  equation: 


s{t)  = 


Cj, 

1 , 


{s,(t)-ER  (t)-ST{t) 


(4) 


where  £•,,  eR,  and  £yare  the  incident,  reflected  and  transmitted  strain  pulses  time  shifted  to  the  front  and  rear  faces 
of  the  specimen,  respectively,  Cb  is  the  sound  speed  in  the  bar  material,  and  /s  is  the  length  of  the  sample.  This 
specimen  strain  rate  is  then  integrated  to  give  the  strain, 


i 

s{t)  =  |  s(t)dt 


(5) 


and  the  transmitted  strain  pulse  is  used  to  calculate  the  reported  one-wave  specimen  stress, 


Ml 

A, « 


(6) 


where  Eb,  and  Ab  are  the  elastic  modulus  and  cross-sectional  area  of  the  bar  material,  respectively,  and  As  is  the 
cross-sectional  area  of  the  sample.  The  two-wave  specimen  stress  is  calculated  using  Equation  6  with  sT 
replaced  by  st  +  sR.  If  true  stress  is  required,  As  is  typically  updated  using  the  strain  calculation,  assuming  that 
volume  is  conserved  during  deformation. 

In  the  quasi-static  experiments,  the  elastic  modulus  was  determined  by  fitting  a  straight  line  to  the  initial,  linear 
part  of  the  stress-strain  curve.  The  yield  stress,  in  both  quasi-static  and  dynamic  experiments,  was  determined  by 
fitting  a  second  order  polynomial  to  the  yield  region  of  the  stress-strain  curve  and  taking  the  derivative  to 
determine  the  maximum. 

RESULTS  AND  DISCUSSION 

Representative  curves  of  the  quasi-static  compression  results  are  shown  in  Figure  2  (a-d)  in  blue.  It  was  found  in 
these  experiments  that  extreme  care  needed  to  be  taken  in  order  to  center  the  samples  on  the  load  path  to  avoid 
shear  failure  in  the  materials.  That  the  measured  stress-strain  curves  in  the  quasi-static  regime  are  serrated  in 
nature  is  believed  to  be  an  artifact  of  the  testing  set-up  rather  than  due  to  any  intrinsic  material  property.  The  tests 
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were  conducted  under  constant  true  strain  conditions  by  continuously  modulating  the  ram  speed  rate  according  to 
the  measured  force  and  ram  displacement,  in  a  feedback  loop.  Factors  such  as  machine  compliance  and  stick- 
slip  loading  on  the  compression  specimen  faces,  along  with  any  visco-elastic  behavior  of  the  epoxy  matrix 
material  can  dramatically  affect  the  control  feedback  to  the  ram,  resulting  in  serrated  loading  curves.  Similar 
serrations  were  not  observed  for  the  SHPB  loading  curves,  which  do  not  use  this  kind  of  feedback  loop  to  effect 
true  strain  rate  conditions.  However,  the  overall  trend  in  the  stress-strain  curve  is  a  rise  to  peak  stress  followed 
by  strain  softening,  more  evident  in  MNML-3  and  5,  which  have  a  lower  volume  fraction  of  particulate,  and 
subsequent  strain  hardening.  This  behavior  is  typical  for  polymer  materials  and  has  been  observed  in  pure  epoxy 
[10].  Additional  testing  using  constant  crosshead  displacement  is  underway  to  confirm  that  the  serrated  nature  of 
the  stress-strain  curve  is  due  to  the  loading  program  rather  than  an  intrinsic  material  behavior. 
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Figure  2:  Stress-strain  curves  for  (a)  MNML-2,  (b)  MNML-3,  (c)  MNML-5,  and  (d)  MNML-8.  The  volume 
average  curves  are  determined  from  a  volume  average  of  the  stress-strain  curves  calculated  from  the 
Mulliken-Boyce  model  fitted  to  epoxy  [10],  the  Zerilli-Armstrong  model  for  99.99%  pure  aluminum  [11]  and 
the  Johnson-Cook  model  for  nickel  200  [12]. 
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Each  material  was  tested  at  2500  Is  using  a  split  Hopkinson  pressure  bar,  with  representative  curves  shown  in 
gold  in  Figure  2  (a-d).  As  discussed  above,  mechanical  equilibrium  in  the  samples  needs  to  be  verified, 
particularly  in  polymer  based  samples.  Figure  3  (a-d)  presents  the  one  wave  and  two  wave  analyses  for 
representative  samples  of  all  four  materials,  as  well  as  the  strain  rate.  It  can  be  seen  from  this  figure  that  all 
materials  are  in  equilibrium  by  5%  strain.  Due  to  the  high  volume  fraction  of  particles,  both  Ni  and  Al,  MNML-2 
shows  a  different  stress-strain  response  than  the  other  composites.  There  is  a  lack  of  strain  softening  post-yield 
and  instead  strain  hardening  is  observed,  which  is  not  present  in  the  other  samples. 

The  yield  stress,  defined  in  the  experimental  set-up  section,  for  each  material  at  both  strain  rates  is  plotted  in 
Figure  4.  If  additional  strain  rates  were  measured,  it  would  most  likely  be  found,  as  is  the  case  with  epoxy,  that 
there  is  a  bilinear  dependance  of  yield  stress  on  strain  rate  due  the  beta  phase  transistion  in  epoxy  [10].  The 
percent  increase  of  the  yield  stress  with  strain  rate  for  the  composites  are  approximately  the  same.  MNML-3 
exhibits  a  larger  increase;  however,  since  there  are  only  two  data  points  its  difficult  to  conclude  if  this  is  real  or 
due  to  variations  in  experimental  measurements.  Additional  testing  at  several  strain  rates  is  underway  to  further 
elucidate  the  strain  rate  dependance  of  yield  stress  in  these  composite  materials. 
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Figure  3:  Comparison  of  one  wave  and  two  wave  stresses  as  well  as  strain  rate  for  (a)  MNML-2,  (b)  MNML- 
3,  (c)  MNML-5  and  (d)  MNML-8. 
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Figure  4:  Strain  rate  dependance  of  yield  stress  for  MNML-2,  3,  5,  and  8. 
Table  2:  Mulliken-Boyce  parameters  for  epoxy  [10] _ _ 


Variable 

Units 

Value 

Ea{e,s) 

MPa 

DMA  data 

Ep{e,e) 

MPa 

DMA  data 

va{0,s) 

0.38 

vp[0,e) 

0.38 

y 0,  a 

s'1 

2.29x1 015 

y  o,/? 

s'1 

2x1 06 

AGa 

J 

3.83x1  O'19 

<1 

J 

3.32x1 0 2° 

&D.a 

0.316 

aP,e 

0.316 

Sss,a 

MPa 

0.58 

ha 

MPa 

300 

CR 

MPa 

14.2 

N 

m"w 

2.3 

A  modeling  approach  was  considered  for  predicting  the  properties  of  the  composites  from  the  properties  of  the 
individual  components.  In  this  method,  a  global  volume  average  for  each  component,  combined  with  a  strength 
model  for  each  component  was  used  to  determine  a  volume  averaged  stress-strain  behavior.  The  volume 
fraction  of  each  component  was  multiplied  by  the  stress  given  by  a  constitutive  equation  for  that  component. 
Each  of  these  was  added  together  to  obtain  a  volume  averaged  result.  For  epoxy,  the  Mulliken-Boyce  model  [13, 
14]  for  glassy  polymers  was  used  with  the  parameters  listed  in  Table  2  [10].  The  equations  defining  this  model 
are  given  in  detail  in  [10,  13]. 

The  Zerilli-Armstrong  model  [11]  was  used  for  pure  aluminum,  with  the  parameters  presented  in  Table  3.  In  this 
model,  based  on  the  thermal  activation  of  dislocation  motion,  the  stress  in  FCC  metals  is  given  by 


<j  =  c0  +  B0£A  exp(-  [3T ) , 


(7) 
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where 


Table  3:  Zerilli-Armstrong  parameters  for  99.99%  pure  aluminum  [11] 


Material 

B0  (MPa) 

Po  ( K'1) 

Pi  HO 

og  (MPa) 

k  (MPa*mm1u) 

Aluminum  (99.99%) 

150 

0.00089 

0.000183 

7 

1.5 

Table  4:  Johnson-Cook  parameters  for  nickel  200  [12] 


Material 

ob  (MPa) 

B  (MPa) 

n 

C 

m 

T melt  (K) 

Nickel  200 

163 

648 

0.33 

0.006 

1.44 

1726 

Cq  =  (7q  +  kl  2 .  (8) 

j3  =  j30-j3Ans,  (9) 

<ris  the  effective  stress,  s  is  the  effective  strain,  where  a  dot  indicates  the  strain  rate,  aG  is  the  stress  contribution 
due  the  initial  dislocation  density,  /  is  the  grain  size,  k  is  the  microstructural  stress  intensity,  and  B0l  Po,  and  Pi  are 
constants.  The  aluminum  particle  size  in  the  composites  (5  or  50  pm)  was  used  for  the  grain  size  term. 

For  the  nickel  component,  the  Johnson-Cook  model  [12]  for  nickel  200  was  used  with  the  parameters  listed  in 
Table  4.  In  this  empirical  model,  multiplicative  terms  are  used  to  describe  the  strain,  strain  rate,  and  temperature 
dependance 


cr  =  cr0  +  Be"  Jl  +  C  In  £■  *][l  -  T  *m 


(10) 


where  a0,  B,  C,  n,  and  m  are  constants, 


with  £q=  1,  and 

T*=  T-Tr  (12) 

Tm  ~Tr 

where  the  subscript  r  refers  to  a  reference  temperature,  typically  room  temperature,  and  the  subscript  m  refers  to 
the  melt  temperature. 

The  stress-strain  curves  calculated  from  volume  averaging  of  the  constitutive  relationships  for  the  component 
materials  are  shown  in  Figure  1,  where  the  low  strain  rate  is  given  by  the  red  curve  and  the  high  strain  rate  is 
given  by  the  green  curve.  For  all  the  composites,  the  comparison  at  low  strain  rates  is  reasonably  good  for  a  first 
approximation.  The  low  strain  rate  behavior  of  MNML-5  is  underpredicted,  which  could  be  due  to  the  equation  of 
state  for  pure  aluminum,  which  was  used.  The  aluminum  powder  used  in  this  study  is  99.7%  pure,  compared  with 
the  99.99%  pure  bulk  aluminum  used  to  develop  the  model.  The  volume  averaged  calculations  underpredict  the 
behavior  in  all  of  the  composites  at  the  high  strain  rates,  indicating  that  there  may  be  interactions  at  the  high  strain 
rates  that  are  not  taken  into  account  purely  by  the  high  strain  rate  behavior  of  the  components. 

The  samples  prepared  in  this  study  were  chosen  as  part  of  a  half-factorial  fractional  design  of  experiments.  As 
such,  analysis  can  be  conducted  to  determine  the  relative  importance  of  each  of  the  three  variables,  aluminum 
particle  size  and  aluminum  and  nickel  volume  fractions,  in  determining  the  mechanical  properties  of  the  composite 
materials.  Table  5  presents  a  summary  of  the  quasi-static  and  dynamic  yield  stresses,  the  elastic  modulus  and 
the  percent  increase  in  stress  from  6x1 0"3  Is  to  2500  Is.  A  comparison  of  the  significance  of  the  three  variables 
tested  on  the  measured  yield  stress,  elastic  modulus  and  strain  rate  slope  revealed  that,  as  far  as  could  be 
determined  by  the  half-factorial  experiments,  they  are  equally  important  in  determining  the  mechanical  properties. 
Typically,  a  half-factorial  design  is  not  attempted  for  a  system  with  only  three  variables.  Additional  work  is 
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underway  to  prepare  samples  for  the  remaining  half  of  the  factorial  design  in  order  to  clarify  the  dependance  of 
the  mechanical  properties  on  these  experimental  variables. 


Table  5:  Summary  of  mechanical  properties  of  composite  materials,  where  Vf  is  the  volume  fraction,  oys  is 
the  yield  strength,  and  E  is  the  elastic  modulus. _ _ _ _ 


Al  Particle  Size 
(pm) 

Vf  Al 

Vf  Ni 

Cys 

(MPa) 

E 

(GPa) 

ays  @  2500/s 
(MPa) 

%  Increase  in 
Stress 

MNML-2 

5 

0.4 

0.1 

97.7 

13.1 

197 

50.4 

MNML-3 

50 

0.2 

0.1 

86.4 

7.6 

197 

56.2 

MNML-5 

50 

0.4 

0 

88.7 

10.3 

182 

51.3 

MNML-8 

5 

0.2 

0 

89.6 

6.0 

188 

52.3 

SUMMARY 

This  paper  presented  the  compressive  stress-strain  behavior  of  epoxy-aluminum-nickel  composites  for  two  strain 
rates,  prepared  using  a  half-factorial  design  of  experiments.  Additional  testing  needs  to  be  conducted  to  totally 
define  the  strain  rate  dependance  of  these  materials.  Two  methods  for  modeling  the  mechanical  response  of  the 
composites  from  the  properties  of  the  constituents  were  used.  In  the  microstructure-based  method,  refined 
models  for  the  component  behavior,  ex.  a  viscoelastic  model  for  epoxy,  need  to  be  included  in  order  to  have 
better  agreement  with  the  experimental  data.  As  a  first  approximation,  the  bulk  volume  averaging  method 
showed  good  agreement  with  the  low  strain  rate  experiments.  However,  there  was  less  agreement  with  high 
strain  rate  experiments  indicating  that  there  may  be  additional  interactions  taking  place  at  high  strain  rates  that 
are  not  accounted  for  by  merely  averaging  the  components.  The  significance  of  the  three  variables  studied  on  the 
mechanical  properties  could  not  be  determined  using  only  a  half-factorial  design  of  experiments  and  additional 
samples  are  being  prepared  to  accomplish  a  full  factorial  design. 
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